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This talk summarizes the recent development in the evaluation of the leading order hadronic contributions to the running of the 
QED fine structure constant a(s), at s = M|, and to the anomalous magnetic moments of the muon (g — 2)^. The accuracy of the 
theoretical prediction of these observables is limited by the uncertainties on the hadronic contributions. Significant improvement has 
been achieved in a series of new analyses which is presented historically in three steps: (I), use of r spectral functions in addition to 
e+e - cross sections, (II), extended use of perturbative QCD and (III), application of QCD sum rule techniques. The most precise 
values obtained are: A« had (M|) = (276.3 ± 1.6) X 10~ 4 , yielding a~ 1 (M^) = 128.933 ± 0.021, and a{f d = (692.4 ± 6.2) X lO" 10 
with which one finds for the complete Standard Model prediction a^ M = (11 659 159.6 ± 6.7) X 10" 10 . For the electron (g - 2) e , the 
hadronic contribution is a£ ad = (187.5 ± 1.8) X 10 -14 . 



1 Introduction 

The running of the QED fine structure constant a(s) 
and the anomalous magnetic moment of the muon are 
famous observables for which the theoretical precision is 
limited by second order loop effects from hadronic vac- 
uum polarization. Both quantities are related via disper- 
sion relations to the hadronic production rate in e + e~ 
annihilation, 



R(s) 



Otot(e + e — > hadrons) 
<T (e+e- — ► 



(1) 



with a (e + e~ — > ^ + aO = 47ra 2 /(3s). While far from 
quark thresholds and at sufficiently high energy y/s, R(s) 
can be predicted by perturbative QCD, theory fails when 
resonances occur, i.e., local quark-hadron duality is bro- 
ken. Fortunately, one can circumvent this drawback by 
using &XeT annihilation data of R(s) and, as proposed 
in Ref. Ill, hadronic r decays benefitting from the largely 
conserved vector current (CVC) , to replace theory in the 
critical energy regions. 

There is a strong interest in the electroweak phe- 
nomenology to reduce the uncertainty on a (M|) which 
used to be a serious limit to progress in the determina- 
tion of the Higgs mass from radiative corrections in the 
Standard Model. Table [l] gives the uncertainties of the 
different Standard Model input expressed as errors on 
sin #wH. Using the former valueB for a(M|), the dom- 
inant uncertainties stem from the experimental sin 2 #w 
determination and from the running fine structure con- 
stant. Thus, any useful experimental amelioration on 
sin 2 #vv requires a better precision of i.e., an im- 

proved determination of its hadronic contribution. 

The anomalous magnetic moment = (g — 2)/2 of 
the muon is experimentally and theoretically known to 
very high accuracy. In addition, the contribution of heav- 
ier objects to a M relative to the anomalous moment of the 



Table 1: Uncertainties of the Standard Model constraints ex- 
pressed in terms of Asin 2 (?w Downward arrows indicate future 
experimental improvement. 



Input 



Asin^vv x 10 5 Uncertainty / Source 



Experiment 


19 


(LEP+SLD)U L 
Aa(M|) = 0J39B 


a(M z ) 


23 


m t 


15 


(CDF+D0)El| Q 


Theory 


5-10 


2-loop EW predictionu 




150 


65 - 1000 GeV 



electron scales as (m M /m e ) 2 ~ 4 x 10 4 . These properties 
allow an extremely sensitive test of the validity of elec- 
troweak theory. The present value from the combined 
H + and jjT measurements □, 



(11 659 230 ±85) x 10" 



(2) 



is expected to be improved to a precision of at n least 
4 x 10~ 10 by the E821 experiment at BraokhavenB (see 
also the contribution to this conference B). Again, the 
precision of the theoretical prediction of a M is limited 
by the contribution from hadronic vacuum polarization 
determined analogously to a(M|) by evaluating a disper- 
sion integral using e + e~ cross sections and perturbative 
QCD. 

2 Running of the QED Fine Structure Constant 

The running of the electromagnetic fine structure 
constant a(s) is governed by the renormalized vac- 
uum polarization function, II 7 (s). For the spin 1 
photon, II 7 (s) is given by the Fourier transform of 
the time-ordered product of the electromagnetic cur- 
rents icm( s ) m the vacuum (q^q v — q 2 g^ v )W 1 (q 2 ) — 
iJ^ze^<0|TO^(s)j^(Q))|0). 
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With Aa(s) = -47raRe [II 7 (s) - EL^O)] and Aa(s) = 
Aaiep(s) + Aahad(s), which subdivides the running con- 
tributions into a leptonic and a hadronic part, one has 



a(s) = 



a(0) 



1 - Aa lcp (s) - Aa had (s) 



(3) 



where 47ra(0) is the square of the electron charge in the 
long- wavelength Thomson limit. 

For the case of interest, s — Mf , the leptonic contri- 

y 



bution at three- loop order has been calculated to bet 

Aai op (M|) = 314.97686 x 1CT 4 . (4) 

Using analyticity and unitarity, the dispersion integral 
for the-contribution from hadronic vacuum polarization 
readslB 



Aa had (M|) 



a(0)M 



3tt 



2- Re / ds 



R(s) 



s(s - M|) - it 



, (5) 



and, employing the identity l/(a:' — x — ie) t ^o = 
P{l/(x'— x)}+iirS(x' — x), the above integral is evaluated 
using the principle value integration technique. 



withx = and/3 M = (l-4m=/s)V2. The 

function if (s) decreases monotonically with increasing s. 
It gives a strong weight to the low energy part of the in- 
tegral (0). About 92% of the total contribution to ajf d 
is accumulated at cm. energies ^/s below 1.8 GeV and 
72% of aj^ ad is covered by the two-pion final state which 
is dominated by the p(770) resonance. Data from vector 
hadmnic r decays published by the ALEPH Collabora- 
tion E3 provide a precise spectrum for the two-pion final 
state as well as new input for the lesser known four-pion 
final states. This new information improves significantly 
the precision of the a„ ad determination!^. 



4 Improvement in 3 Steps 

A very detailed and rigorous evaluation of both a(M|) 
and (g — 2) M was performed by S. Eidelman and 
F. Jegerlehner in 1995 H which since then is frequently 
used as standard reference. In their numerical calcu- 
lation of the integrals (||) and ([?]), the authors use ex- 
clusive e + e~ — > hadrons cross section measurements be- 
low 2 GeV cm. energy, inclusive R measurements up to 
40 GeV and finally perturbative QCD for above 40 GeV. 
Their results to which I later will refer are 



3 Muon Magnetic Anomaly 

It is convenient to separate the Standard Model predic- 
tion for the anomalous magnetic moment of the muon, 



(g — 2)u/2, into its different contributions, 



fl SM = a QED 



,had 



,wcak 



(6) 



where a® ED = (11658470.6 ± 0.2)^x 10" 10 is the 
pure electromagnetic contribution (seeL3 and references 
therein), a^ ad is the contribution from hadronic vacuum 



polarization, and a™ ak = (15.1 ± 0.4) x 10~ 10 fc3 
accounts for corrections due to exchange of the weak in- 
teracting bosons up to two loops. 

Equivalently to Aahad(-^z), by virtue of the analyt- 
icity of the vacuum polarization correlator, the contribu- 
tion of the hadronic vacuum polarization to a M can be 
calculated via the dispersion integralEj 



had _ " 2 (0) f ,K{ S ) 



ds R(s) , 



where K(s) denotes the QED kernelE] , 



(7) 



K(s) =x i [l-—J+ (1 + xY + - 

x2 \ , ( 1 + 3J )„2 



x ln(l + x) + 



2 J (1-x) 



x z lnx , (8) 



Aa had (M|) = (279. 7± 6.5) x 10" z 
= (702.4 ± 15.3) x 10" 



^had 



(9) 



Due to improvements on the experimental sides these 
theoretical evaluations are insufficient for present needs. 
Fortunately, new data and a better understanding of the 
underlying QCD phenomena led to new and significantly 
more accurate determinations of the hadronic contribu- 
tions to both observables. 

(I) Addition of Precise r Data 

Using the conserved vector current (CVC) it was shown 
in Ref. EJ that the addition of precise r spectral func- 
tions, in particular of the r~ — ► tt^tt v t channel, to the 
e + e~ annihilation cross section measurements improves 
the low-energy evaluation of the integrals (||) and (0). 
Hadronic r decays into ud' isovector final states occur 
via exchange of a virtual W~ boson and have therefore 
contributions from vector and axial- vector currents. On 
the contrary, final states produced via photon exchange 
in e + e~ annihilation are always vector but have isovec- 
tor and isoscalar parts. The CVC relation between the 
vector two-pion r spectral function Wj=i(t — * irir v T ) 
and the corresponding isovector e + e" cross section at 
energy-squared s reads 



(e+e- 



+ 47ra 2 (0) o 

7T 7T ) = Vj = l(T — » 7T7T V T ) 
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where vj—i(t — > 7r7r°j/ r ) is essentially the hadronic 
invariant mass spectrum normalized to the two-pion 
branching ratio and corrected by a kinematic factor ap- 
propriate to t decays with spin J = 1 EZL The two-pion 
cross sections (incl. the r contribution) in different en- 
ergy regions are depicted in Fig. [l[ Excellent agreement 
between r and e + e~ data is observed. For the four pion 
final states, isospin rotations must be performed to re- 
late the respective r charges to the corresponding e + e~ 
topologiesO. 

Effects from SU (2) violation 

Hadronic spectral functions from r decays are directly 
related to the isovector vacuum polarization currents 
when isospin invariance (CVC) and unitarity hold. For 
this purpose one has to worry whether the breakdown 
of CVC due to quark mass effects (m u ^ md generat- 
ing <9 M J M ~ (m u — rrid) for a charge-changing hadronic 
current J M between u and d quarks) or unknown isospin- 
violating electromagnetic decays have non-negligible con- 
tributions within the present accuracy. Expected de- 
viations from CVC due to so-called second class cur- 
rents as, e.g., the decay r~ — * ix~r]v T where the cor- 
responding e + e~ final state Tr°n (C=+l) is strictly for- 
bidden, have estimated branching fractions of the order 
of (m u — md) 2 — E3, while the experimental uit 

per limit amounts to B(r — > ix~r)v T ) < 1.4 x 10 — 4 113 
with 95% CL. SU (2) symmetry breaking caused by elec- 
tromagnetic interactions can occur in the p ± -p° masses 
and widths. Hadronic contributions to the p ± -p° width 
difference are expected to be much smaller since they 
are proportional to (m u — md) 2 . The total expected 
SU(2) violation in the p width is estimated in Ref. tl to 
be (r p ± - I»/T p = (2.8 ± 3.9) x l(r 3 , yielding the 
corrections 



Sa had = -(1.3 ±2.0) x 10 



-10 



SAa^M'i) = -(0.09 ± 0.12) x 10~ 4 , (10) 

to the respective dispersion integrals when using the 
r~ — > 7r _ 7r°^ T spectral function in addition to e + e~ data. 

Evaluation of the dispersion integrals (j^) and (Q) 

Details about the non-trivial task of evaluating in a co- 
herent way numerical integrals over data points which 
have statistical and correlated systematic errors between 
measurements and experiments are given in Ref. □. The 
procedure is based on an analytical x 2 minimization, tak- 
ing into account all initial correlations, and it provides 
the averages and the covariances of the cross sections 
from different experiments contributing to a certain final 
state in a given range of cm. energies. One then ap- 
plies the trapezoidal rule for the numerical integration of 
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Figure 1: Two-pion cross section as a function of the cm. energy- 
squared. The band represents the result of the averaging procedure 
described in the text within the diagonal errors. The lower left hand 
plot shows the chiral expansion of the two-pion cross section used 
(see Ref.lil). 



the dispersion integrals (||) and (Q), i.e., the integration 
range is subdivided into sufficiently small energy steps 
and for each of these steps the corresponding covariances 
(where additional correlations induced by the trapezoidal 
rule have to be taken into account) are calculated. This 
procedure yields error envelopes between adjacent mea- 
surements as depicted by the shaded bands in Fig. 

Results 

With the inclusion of the r vector spectral functions, the 
hadronic contributions to the running a(M|) and to a p 
are found to beB 



Aa had (M z 2 ) = (281.0 ±6.2) x KT 4 , 
„had _ (jQl.l ± 9.4) x 10~ 10 



(11) 



with an improvement for a p of about 40% compared to 
the previous evaluation (|^) , while there is only a marginal 
improvement of Aa (M|) for which the dominant uncer- 
tainties stem from higher energies. 

(II) Extended Theoretical Approach 

The above analysis shows that in order to improve the 
precision on Aa(Af|), a more accurate determination of 
the hadronic cross section between 2 GeV and 10 GeV 
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is needed. On the experimental side there are ongo- 
ing R jaeasurements performed by the BES Collabo- 
ration E3. On the other hand, QCD ^nalyses using 
hadronic jE_dacays performed by ALEPHl^I and recently 
by OPALcjO revealed excellent applicability of the Wil- 
son Operator Eroduct Expansion (OPE) CJ (also called 
SVZ approach Ej), organizing perturbative and nonper- 
turbative contributions to a physical observable using the 
concept of global quark-hadron duality, at the scale of the 
r mass, M T ~ 1.8 GeV. Using moments of spectral func- 
tions, dimensional nonperturbative operators contribut- 
ing to the t hadronic width have been fitted simultane- 
ously and turned, out to be small. This encouraged the 
authors of Ref. c3 to apply a similar approach based on 
spectral moments to determine the size of the nonpertur- 
bative contributions to integrals over total cross sections 
in e + e~ annihilation, and to figure out whether or not 
the OPE, i.e., global duality is a valid approach at rela- 
tively low energies. 

Theoretical prediction of R(s) 

The optical theorem relates the total hadronic cross sec- 
tion in e + e~ annihilation, R(sq), at a given energy- 
squared, so j to the absorptive part of the photon vacuum 
polarization correlator 



R(s ) = 127rimn(s + ie) 



(12) 



Perturbative QCD predictions up to next-to-next-to 
leading order as well as second order quark mass cor- 
rections far from the production threshold and the first 
order dimension D — 4 nonperturbative term are avail- 
able for the Adler D-functionEj, which is the logarithmic 
derivative of the correlator II, carrying all physical infor- 
mation: 



D(s) = -12tt 2 s 



dIL(s) 
ds 



This yields the relation 



R( So ) = J_ <f ±D( S ) 
2m / s 



(13) 



(14) 



|s|=s 



where the contour integral runs counter-clockwise around 
the circle from s =_S(>-Jpje to s = sq — ie. The Adler 
function is given byEa'Ej'ES 



D fi (-s)=N c J2Q 2 f{ 



i+do^i+d, f^i) 2 +d 2 (^: 3 

IT \ 7T / \ IT 



6 + 28- 



(295.1 - 12. 3 n/) 



2tt 2 / Ua a (s)\ (%-GG) 

3 v is 7T y S 2 

8tt 2 (l - ( m /g/g/) 




(15) 



where additional logarithms occur when /i 2 ^ s and 
H being the renormalization scale^]. The coefficients 
of the massless perturbative part are do = 1, di = 
1.9857 - 0.1153 n/, d 2 = d 2 + /3 2 tt 2 /48 with fa = 
11 - 2n//3, d 2 = -6.6368 - 1.2001 n f - 0.0052 n) - 
1.2395 {J2 f Qf) 2 /NcJ2f Q) and n f being the number 
of involved quark flavours. The nonperturbative opera- 
tors in Eq. (l||) are the gluon condensate, {(a s /ir)GG) , 
and the quark condensates, (m/q/q/). The latter obey 
approximately the PCAC relations 



(m« + m d )(uu + dd) ~ -2/ 2 m 2 , 

m s (ss) ~-/ 2 (m^-m 2 ) , (16) 

with the pion decay constant f„ = (92.4 ± 0.26) MeV@. 
In the chiral limit the equations / w = Jk and (uu) — 
(dd) = (ss) hold. The complete dimension D = 6 and 
D = 8 operators are parametrized phenomenologically in 
Eq. ( |l5|) using the saturated vacuum expectation values 
(Oe) and (Og), respectively. 

Although the theoretical prediction of R using 
Eqs. (14) and ( |l5|) assumes local duality and therefore 
suffers from unpredicted low-energy resonance oscilla- 
tions, the following integration, Eqs. (|))/(0), turns du- 
ality globally, i.e., the nonperturbative oscillations are 
averaged over the energy spectrum. However, a system- 
atic uncertainty is introduced through the cut at explic- 
itly 1.8 GeV so that non- vanishing oscillations could give 
rise to a bias after integration. The associated systematic 
error is estimated in Ref. Ell by means of fitting different 
oscillating curves to the data around the cut region, yield- 
ing the error estimates A(Aahad(-^D) = 0-15 x 10~ 4 
and AaJ^ ad = 0.24 x 10~ 10 , from the comparison of the 
integral over the oscillating simulated data to the OPE 
prediction. These numbers are added as systematic un- 
certainties to the corresponding low-energy integrals. 

In. asymptotic energy regions we use the formulae of 
Ref. which include complete quark mass corrections 
up to order a 2 to evaluate the perturbative prediction of 
R(s) entering into the integrals (||) and (0). 

a The negative energy-squared in D(-s) of Eq. @ is introduced 
when continuing the Adler function from the spacclike Euclidean 
space, where it is originally defined, to the timelike Minkowski 
space by virtue of its analyticity property. 
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Theoretical uncertainties 

Details about the parameter errors used to estimate the 
uncertainties accompanying the theoretical analysis are 
given in Refs. c3S3. Theoretical uncertainties arise from 
essentially three sources 

(i) The perturbative prediction. The estimation of the- 
oretical errors of the perturbative series is strongly 
linked to its truncation at finite order in a s . 
This introduces a non-vanishing dependence on 
the choice of the renormalization scheme and the 
renormalization scale. Furthermore, one has to 
worry whether the missing four-loop order contri- 
bution dz{a s /n) A: gives rise to large corrections to 
the perturbative series. An additional uncertainty 
stems from the ambiguity between the results on R 
obtained using contour-improved fixed-order per- 
turbation theory (FOPT C i) and FOPT only (see 
Ref.B). 

(ii) The quark mass correction. Since a theoretical 
evaluation of the integrals @ and (Q) is only ap- 
plied far from quark production thresholds, quark 
mass corrections and the corresponding errors are 
small. 



(Hi) 



The nonperturbative contribution. In order to de- 
tach the measurement from theoretical constraints 
on the nonperturbative parameters of the OPE, 
the dominant dimension D = 4, 6, 8 terms are de- 
termined experimentally by means of a simultane- 
ous fit of weighted integrals over the inclusive low 
energy e + e~ cross section, so-called spectral mo- 
ments, to the theoretical prediction obtained from 
Eq. (|l5|). Small uncertainties are introduced from 



possible deviations from the PC AC relations (16) 



The spectral moment fit of the nonperturbative op- 
erators results in a very small contribution from the 
OPE power terms to the lowest moment at the scale of 
1.8 GeV (repeated and ccpficmed at 2.1 GeV), as ex- 
pected from the r analysesEEl The value of (^GG) = 
(0.037 ± 0.019) GeV 4 found for thc^lJuon condensate 
is compatible with other evaluations Ejcj. The analysis 
proved that global duality holds at 1.8 GeV and nonper- 
turbative effects contribute only negligibly, so that above 
this energy perturbative QCD can replace the rather im- 
precise data in the dispersion integrals (||) and (^). 

Results 

The R(s) measurements and the corresponding theoret- 
ical prediction are shown in Fig. |^. The shaded bands 
indicate the regions where data are used instead of the- 
ory to evaluate the dispersion integrals. This is below 




Figure 2: Inclusive hadronic cross section ratio in e+e - anni- 
hilation versus the cm. energy y^s. Additionally shown is the 
QCD prediction of the continuum contribution as explained in the 
text. The shaded areas indicate regions were experimental data 
are used for the evaluation of Acn la ^(M^) and ajj ad in addition to 
the measured narrow resonance parameters. The exclusive e + e — 
cross section measurements at low cm. energies are taken from 
DM1,DM2,M2N,M3N,0LYA,CMD,ND and r data from ALEPH 
(sec Ref. EJ for detailed information). 



1.8 GeV and at cc threshold energies. Good agreement 
between data and QCD is found above 8 GeV, while at 
lower energies systematic deviations are observed. The R 
measurements in this region are essentially provided by 
the 772E3 and MARK I\B collaborations. MARK I data 
above 5 GeV lie systematically .above the measurements 
of the Crystal BallE3 and MDlN Collaborations as well 
as above the QCD prediction. 

The combination of the theoretical and experimental 
evaluations of the integrals yields the resultsc3 



Aa had (Mi) 

„had 



(277.8 ±2.2 oxp ±1.4 thco ) x 10~ 4 , 
(695.1 ±7.5 oxp ±0.7 t hco) x KT 10 , (17) 



with a significant improvement by more than a factor 
of two for Aa(Af|), and a 20% better accuracy on a^ ad 
compared to the numbers iO). 

The authors of Ref. EBimproved the above anal- 
ysis in the charm region by normalizing experimen- 
tal results in the theoretically not accessible region (at 
least locally) so that they match perturbative QCD at 
safe energies below and above the occurrence of reso- 
nances. The so-renormalized data show excellent agree- 
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ment among different experiments which supports the 
hypothesis made that experimental systematic errors are 
completely correlated over the whole involved energy 
regime. The result (after_correcting for the small top 
quark contribution) reads Ej 



Aa had (M|) = (276.7 ±1.7) x 10" 



(18) 



Another, very elegant method based on an analytical 
calculation of the unsubtracted dispersion relation, corre- 
sponding to the subtracted integral (H), was presented in 
Rcf. E3. Only the low-energy pole contribution is taken 
from data, while the contribution from higher energies 
is calculated analytically using the two-point correlation 
function given in Ref. E 2 ], and the renormalization group 
equations forJihe running quantities. This leads to the 
precise result c3 



Aa had (M|) = (277.2 ±1.9) x 10" 



(19) 



Both numbers 
Aa(M|) from Eq. 




are in agreement with 



(III) Constraints from QCD Sum Rules 

It was shown in Ref. that the the previous determi- 
nations can be further improved by using finite-energy 
QCD sum rule techniques in order to access theoret- 
ically energy regions where locally perturbative QCD 
fails. This idea was first presented in Ref. ell. In principle, 
the method uses no additional assumptions beyond those 
applied in Section ^j. However, parts of the dispersion 
integrals evaluated at low-energy and the cc threshold 
are obtained from values of the Adler D-function itself, 
for which local quark-hadron duality is assumed to hold. 
One therefore must perform an evaluation at rather high 
energies (3 GeV for u, d, s quarks and 15 GejLfor the c 
quark contribution have been chosen in Ref. EU) to sup- 
press deviations from local duality due to nonperturba- 
tive phenomena. 

The idea of the approach is to reduce the data con- 
tribution to the dispersion integrals by subtracting ana- 
lytical functions from the singular integration kernels in 
Eqs. (0) and (Jfy, and adding the subtracted part sub- 
sequently by using .theory only. Two approaches have 
been applied in Ref.Ei: first, a method based on spectral 
moments is defined by the identity 



F= I dsR(s)[f(s)-p n (s)} 
2iri J s 

|s|=80 



Pn(so)-Pn(s)]D uds (s) , (20) 



with P n (s) = f*dtp n (t) and f(s) = a(0) 2 K {s)/(3n 2 s) 



for F = a^ 2mKj v/53] , as well as f(s) = a(0)M 2 /(3ns(s- 
A/f)) for F = Aa had (M|) [2rn7ri The regular func- 

tions p n {s) approximate the kernel /(s) in order to re- 
duce the contribution of the first integral in Eq. (20) 
which has a singularity at s — and is thus evaluated 



using experimental data. The second integral in Eq. (20) 



can be calculated theoretically in the framework of the 
OPE. The functions p n (s) are chosen in order to reduce 
the uncertainty of the data integral. This approximately 
coincides with a low residual value of the integral, i.e., a 
good approximation of the integration kernel f(s) by the 
p n {s) defined as Ell 



(21) 



with the form (1 — s/sq) in order to ensure a vanishing 
integrand at the crossing of the positpe real axis where 
the validity of the OPE is questionedc3. Polynomials of 
order s n involve leading order nonperturbative contribu- 
tions of dimension D — 2(n+l). The analysis is therefore 
restricted to the linear n = 1 case only. 

A second approach uses the dispersion relation of the 
Adler D-function 




oo 

Df{Q 2 ) = Q 2 J ds 



(s + Q 2 ) 2 



(22) 



for space- like Q 2 = —q 2 and the quark flavour /. The 
above integrand approximate the integration kernels in 
Eqs. (||) and (0), so that the modified Eq. ( p0[ ) reads 



F = J ds R 

4m 



Data 



/(*)- 



+ A F [D uds {Q 2 )-Q 2 Jds 



A F Q 2 ' 
(.s + Q 2 ) 2 . 



Q 2 



(23) 



with a normalization constant Ap to be optimized for 
both Aa(M|) and a£ ad . 

The theoretical errors are evaluated correspondingly 
to the analysis presented in Section |[ 

Results 

A \ 2 fit taking into account the experimental and the- 
oretical correlations between the polynomial moments 
yields for the first (spectral moment i-approach (hadronic 
contribution from 2m w to 1.8 GeV) Ell 



Aa had (M|) [1 . 8 GcV ] 



had 
a M, [1.8 GeV] 



(56.53 ± 0.73ex P ± 0.39 th co) x 10" 
(634.3 ± 5.6 exp ± 2.1 t heo) x 1CT 10 
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Table 2: Contributions to Aa^ sl ^(M^), aj^ ad and to aj? ad from the different energy regions. The subscripts in the first column give the 
quark flavours involved in the calculation. 



Energy (GeV) 


Aahad (M|) x 10 4 


a h ad x 1Q io 


oT d x io 14 


(2m T - l.8)uds 


56.36 ± 0.70 cxp ±0.18thco 


634.3 ± 5.6 exp ± 2.1 th Bc 


173.67 ± 1.7exp ± 0.6theo 


(1.8 - 3.700) uda 


24.53 ± 0.28 thc o 


33.87 ± 0.46thco 


8.13 ± O.lltheo 


V>(1S* 2S,3770) C + (3.7-5)„ dsc 


24.75 ± 0.84 cxp ± 0.50thco 


14.31 ± 0.50 cxp ± 0.21 thc o 


3.41 ±0.12 cxp ±0.05 th oo 


(5 - 9.3) u( j S c 


34.95 ± 0.29 thc o 


6.87 ± O.lltheo 


1.62 ± 0.03 thco 


(9-3 ^-2)udscb 


15.70 ± 0.28 thc o 


1.21 ± 0.05 tho o 


0.28 ± 0.02 th eo 


(12 - Oo)udscb 


120.68 ± 0.25 thc o 


1.80 ± O.Olthoo 


0.42 ± 0.01 t h C o 


(2m t - oo) t 


-0.69 ± 0.06 th oo 


« 


« 


(2m T - oo) u dscbt 


276.3 ± l.lexp ± l.lthco 


692.4 ± 5.6 cxp ± 2.6th co 


187.5 ± 1.7 oxp ± 0.7 t heo 



while the. dispersion relation approach gives {y/Q^ = 
3 GeV) El 



AQ h ad(A/|; 



[1.8 GeV] 



had 
a M, [1.8 GeV] 



(56.36 ± 0.70ex P ± 0.18 t heo) X 10" 
(632.5 ± 6.2 cxp ± 1.6 t heo) X 10" 10 



Only the most precise of the above numbers are used 
for the final results. The above theory-improved results 
can be compared to the corresponding pure experimen- 
tal values, Aa had (M|) [1 . 8 GcV ] = (56.77 ± 1.06) x 1(T 4 
and a^f l s GcV] = (635.1 ± 7.4) x 1CT 10 , showing clear 
improvement. 

For the charm threshold region only the ap- 
proach ( p3| ) is used giving for the hadronic_contributions 
from 3.7 GeV to 5 GeV {yfCf = 15 GeV)0 



Aa?had(A^z)[3.7, 5 GeV] 
had 

a fi, [3.7, 5 GeV] 



(24.75 ± 0.84 cxp ± 0.50thco) x 10" 4 , 
(14.31 ± 0.50 cxp ± 0.21 t hco) x TO" 10 



for which compared with the pure data results, 
Aa had (Ml) [3.7, 5 GcV] = (25.04 ± 1.21) x 10~ 4 and 
(14.44 ± 0.62) x lO" 10 , only a slight im- 



7 had 

>, [3.7, 5 GeV] 



provement is observed. 
5 Final Results 

Table || shows the experimental and theoretical evalu- 
ations of Aahad(M|), a^ ad and a\ a,A for the respective 
energy regimes^. Experimental errors between different 
lines are assumed to be uncorrelated, whereas theoreti- 
cal errors, but those from cc and bb thresholds which are 
quark mass dominated, are added lincarily. 

According to Table the combination of the theo- 
retical and experimental evaluations of the integrals @ 
and (Q) yields the final results 

Aa had (Af z 2 ) = (276.3 ± l.lexp ± l.l theo ) x 10~ 4 , 
a- 1 (Ml) = 128.933 ±0.015 exp ±0.015 theo , 

b The evaluation of a had follows the same procedure as ajj ad . 



,had 



,SM 



(692.4 ±5.6 cxp ±2.6 the o) X 10" 



o-- = (11 659 159.6 ± 5.6 cxp ± 3.7 th co) x 10 



-10 



and af d = (187.5±1.7 cxp ±0.7 t hco) X 10~ i4 for the lead- 
ing order hadronic contribution to a e . The improvement 
for a(M|) compared to the previous results (|T?f) amounts 



to 40% and a^ ad is about 17% more precise than (17) 



The total of. value includes an additional contri- 
bution from non-leading ordex-hadronic vacuum polar- 
ization summarized in Refs. Ou to be et^ ad [(a/7r) 3 ] = 
(-10.0 ± 0.6) x 10~ 10 . Also the light-by-light scatter- 
ing (LBLS) contribution has recently been reevaluated 



m Refs.l 4 ^ andcJ of which the average (aJ} ad [LBLS]) = 
(-8.5 ± 2.5) x 10~ 10 is used here. 

Figures || and |] show a compilation of published re- 
sults for the hadronic contributions to a(M|) and a M . 
Some authors give the contribution for the five light 
quarks only and add the top quark part separately. This 
has been corrected for in Fig. ||. 

The new result for a(s) is exploited in Ref. Ej to 
repeat the global electroweak fit in order to adjust the 
mass of the Standard Model Higgs boson, MHi ggs . The 
prediction of the Standard Model is obtained from the 
ZFITTER electroweak libraryEJ. The fit yields M H ig g8 = 



83 
CL 



eV and an upper mass limit of 202 GeV at 95% 



6 Conclusions 

This note summarizes the recent effort that has been 
undertaken in order to ameliorate the theoretical predic- 
tions for a(Mj) and a^ ad , crucially necessary to main- 
tain the sensitivity of the diverse experimental improve- 
ments on the Standard Model Higgs mass, on the one 
hand, and tests of the electroweak theory on the other 
hand. The new value of a' 1 {Ml) = 128.933 ± 0.021 
for the running fine structure constant proves as be- 
ing sufficiently accurate to unaffect the global Standard 
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Model fit. On the contrary, more effort is needed to 
further improve the precision of the hadronic contribu- 
tion to the anomalous magnetic moment of the muon 
below the intended experimental accuracy at BNL-E821, 
which is about 4 x 1CP 10 0. Fortunately, new low en- 
ergy data are expected in the near future from r decays 
(CLEO, OPAL, DELPHI, BABAR) and from e+e~ an- 
nihilation (BES, CMD II, DA<FNE). Additional support 
might come from the theoretical side using chiral pertur- 
bation theory to access the low energy inverse moment 
sum rules (^|) and (Q) in applying^ e.g., a similar proce- 
dure as the one suggested in Ref.El 
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